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ABSTRACT 

Aims. In this work we re-investigated the empirical relation between the hard X-ray peak flux and the outburst waiting time found 
previously in the black hole transient GX 339-4. We tested the relation using the observed hard X-ray peak flux of the 2007 outburst 
of GX 339-4, clarified issues about faint flares, and estimated the lower limit of hard X-ray peak flux for the next outburst. 
Methods. We included Swift/BAT data obtained in the past four years. Together with the CGRO/BATSE and RXTE/HEXTE light 
curves, the observations used in this work cover a period of 18 years. 

Results. The observation of the 2007 outburst confirms the empirical relation discovered before. This strengthens the apparent link 
between the mass in the accretion disk and the peak luminosity of the brightest hard state that the black hole transient can reach. We 
also show that faint flares with peak fluxes smaller than about 0.12 crab do not affect the empirical relation. We predict that the hard 
X-ray peak flux of the next outburst should be larger than 0.65 crab, which will make it at least the second brightest in the hard X-ray 
since 1991. 
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1. INTRODUCTION 



GX 339-4 is a black hole transient discovered more than 30 
years ago. It has a mass function of 5.8 M R , a low mas s com- 
panion star and a distance of > 7 kpc dMarkert et all 1973t 
iHvnes et al.L [2003 1 Sha hbaz et all l2001b IZdziarski et aTfl2004l) 
It is o ne of the black hole transients with the most frequent out- 
bursts dKong et al.Ll2002tlZdziarski et al.Ll2004T) . lYu et al.l(l2007l) 
analyzed the long-term observations of GX 339-4 made by the 
Burst and Transient Source Experiment (BATSE) on board the 
Compton Gamma-Ray Observatory (CGRO) and the Rossi X- 
ray Timing Explorer (RXTE) since May 31, 1991 until May 23, 
2005. They found a nearly linear relation between the peak flux 
of the low/hard (LH) spectral state that occurs at the beginning 
of an outburst and the outburst waiting time defined based on 
the hard X-ray flux peaks. The empirical relation indicates a link 
between the brightest LH state that the source can reach and the 
mass stored in the accretion disk before an outburst starts. 

After then the source underwent an outburst in 2007. The 
2007 outburst and any future outbursts can be used to test and 
refine the empirical relation. Here we show that the hard X-ray 
peak flux of the 2007 outburst falls right on the empirical relation 
obtained bv I Yu et all d2007l) . proving that the empirical relation 
indeed holds. By including the most recent monitoring observa- 
tions with the Swift/BAT in the past four years, we re-examine 
the empirical relation and make a prediction for the hard X-ray 
peak flux of the next bright outburst for a given waiting time. 
We also clarify issues related to faint flares that have been seen 
in the recent past. 



2. OBSERVATION AND DATA ANALYSIS 

We made use of observations performed with BATSE (20- 
160 keV) covering from May 31, 1991 to May 25, 2000, 
HEXTE (20-2 50 keV) c o vering from January 6, 1996 to January 
2, 2006, as in lYu et al.l (120071 and recent monitoring results 
of Swift/BAT that are publicly available (15-50 keV) cover- 
ing from February 13, 2005 to August 31, 2009. The BATSE 
data were obtained in crab unit. The fluxes of the Crab were 
305 counts s _1 and 0.228 counts s" 1 cnr 2 for HEXTE and BAT 
respectively. These values were used to convert the sou rce fluxes 
into t he unit of crab. Following the previous study dYu et all 
2007), the light curves were rebinned to a time resolution of 
10 days. It is worth noting that the X-ray fluxes quoted below 
all correspond to 10-day averages, including those obtained in 
the empirical relation and the predicted fluxes. 

The combined BATSE, HEXTE and BAT light curves are 
shown in FigQ] The triangles marked with 1-8 indicate the ini- 
tial hard X-ray peaks during the rising phases of the outburst 
1-8, and those with 5 e -8 e indicate the ending hard X-ray peaks 
during the decay phases of the outburst 5-8. Outburst 1-7 were 
studied in lYu et all (120071) . Outburst 8 is the 2007 outburst that 
occurred after the empirical relation was obtained. The waiting 
time of outburst 8 is determined in the same way as in the pre- 
vious study, i.e., the time separation between the peaks 7 e and 
8 and the peak 7 e is the hard X-ray peak associated with the 
HS-to-LH transition. In order to show how the peaks are cho- 
sen, we also plotted the soft X-ray light curves obtained with the 
RXTE/ASM and the hardness ratios between the ASM and the 
BATSE or HEXTE or BAT fluxes in Fig [2] This explicitly shows 



2 



Wu et al.: On the Relation of Hard X-ray Peak Flux and Waiting Time in GX 339-4 (RN) 



1.2 

1.0 

S 0.8 

~b 0.6 
fa" 0.4 

0.2 

0.0 



6-5 e< 



,^8-7„ 



BATSE (20-160 keV)- 
HEXTE (20-250 kcV 
BAT (15-50 kcV) - 



^5-4 



200 400 600 800 1000 1200 
Waiting Time (day) 




1000 2000 3000 4000 5000 6000 7000 

TJD-8000 (Day) 

Fig. 1. The long-term hard X-ray light curves of GX 339- 
4 from the observations with CGRO/BATSE (empty circles), 
RXTE/HEXTE (filled circles) and Swift/BAT (squares). The in- 
set panel is the relation between the LH state peak flux and the 
waiting time following Yu et al. (2007). The dashed-dotted line 
indicates the flux level of 0.1 crab, under which X-ray flares 
appears not to affect the empirical relation. The triangles indi- 
cate the LH state peaks used to calculate the waiting times. The 
dashed line in the inset panel shows the best-fitting linear model. 



that the hard X-ray peaks at the end of outbursts correspond to 
the HS-to-LH state transitions. The initial hard X-ray peak, on 
the other hand, is normally the first prominent one during the 
initial LH state. Due to the hys teresis effect of spectral state tran- 
sitions (Miyam oto et all fl 9951) . the source would have very low 
luminosity after the HS-to-LH transition during the outburst de- 
cay. We took the hard X-ray peak corresponding to the HS-to-LH 
state transition such as peak 7 e , as the end of the previous out- 
burst, i.e., the starting time to calculate the waiting ti me of the 
follow ing outburst (see the definition of waiting time in lYu et all 
120071) . 

Due to the relatively low sensitivity of BATSE, flares with 
10-day averaged peak flux at or below about 0.1 crab could not 
be identified as individual outbursts. It is therefore worth noting 
that the current empirical relation is determined based on out- 
bursts with hard X-ray peak fluxes above about 0.2 crab. In re- 
cent years with more sensitive observations of Swift/BAT, we 
have observed several faint flares in this source. These flares 
would not have been clearly seen in the BATSE 10-day averaged 
light curve and would not have been taken as single outbursts if 
BATSE had operated. Therefore we ignored these flares although 
they were clearly seen with Swift/BAT. We will discuss the faint 
flares later on. 

We found that the data point of outbu rst 8 follows the em- 
pirical relation reported in lYu et aU (120071) . as shown in the inset 
panel of FigQ] The deviation from the empirical relation is only 
-0.034 crab. The linear Pearson's correlation coefficient for all 
the 7 data points is 0.997, again indicating a nearly linear rela- 
tion between the hard X-ray peak flux F p and the waiting time 
T w . 

A linear fit to this relation gives F p = (9.25+0.06)xlO~ 4 T w - 
(0.039 + 0.005), where F p is in unit of crab and T w in units 
of days. This updated rel ation is almost identical to the one re- 
ported in I Yu et aU 12007). The intrinsic scattering of the data is 
0.014 crab, which defines a +0.014 crab bound of the linear rela- 
tion. The intercept of the best-fitting linear model on the waiting 
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Fig. 2. Top: the long-term hard X-ray light curves of GX 339- 
4 from the observations with CGRO/BATSE (empty circles), 
RXTE/HEXTE (filled circles) and Swift/BAT (squares). Middle: 
the long-term soft X-ray fight curve of the RXTE/ASM (2- 
10 keV). Bottom: The hardness ratios between the hard X- 
ray fluxes and the ASM fluxes: BATSE/ASM (empty circles), 
HEXTE/ASM (filled circles) and BAT/ASM (squares). The ini- 
tial hard X-ray peaks and the hard X-ray peaks at the end of the 
HS-to-LH state transition in the outburst decays are marked with 
dotted lines and solid lines, respectively. The HS and LH state 
hardness ratio thresholds were set to 4 and 0. 1 respectively, as 
indicated by the dashed lines. 



time axis is T w = 42 days when F p = crab. Considering the 
intrinsic scattering and the model uncertainty, we obtained an 
intercept T w = 42 + 20 days. This means that the hard X-ray 
peak of any outburst should be at least 42 + 20 days after the end 
of the previous outburst, which is determined as the hard X-ray 
peak corresponding to the HS-to-LH transition. 

The refined empirical relation enables us to approximately 
estimate the hard X-ray peak flux (10-day average) for the next 
bright outburst in GX 339-4. The updated relation gives the peak 
flux of the next bright outburst as F pn = 9.25 x 10~ 4 (Dayo9 + 
Trfse) + 0.44 crab, where Dayo9 is the number of days in 2009 
when a future outburst starts and T rise is the rise time in unit of 
day for the next outburst to reach its initial hard X-ray peak. The 
hard X-ray peak flux can be predicted almost as soon as the next 
outburst occurs because the rise time is nearly a small constant 
compared with the waiting time. The source has remained inac- 
tive for about 750 days since the end of the 2007 outburst. This 
gives that the hard X-ray peak flux of the next outburst should 
be at least 0.65 crab (Fig [3}, making it the second brightest out- 
burst since 1991, brighter than all the outbursts except outburst 
6. Again notice that only for an outburst brighter than about 
0.12 crab can such a prediction be made based on the empir- 
ical relation. We have shown that the empirical relation holds 
if faint hard X-ray flares are ignored. For example the flare of 
about 0.08 crab in March 2006 does not affect the peak flux of 
the 2007 outburst. This suggests that the flare of about 0. 1 crab 
in March 2009 will not affect the hard X-ray peak flux of next 
bright outburst significantly. 

The negligible effect of the faint flares on the empirical rela- 
tion is also consistent with the consideration of the actual value 
range of T w . The intersection of the best-fitting linear empirical 
relation on the time axis indicates that the hard X-ray peak of a 
major, bright outburst must occur more than 42 + 20 days after 
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Fig. 3. Prediction of the hard X-ray peak flux for next bright out- 
burst based on the updated empirical relation. The dashed line 
indicates the regime excluded for the next bright outburst. The 
arrow indicates the lower limit of the hard X-ray flux of the next 
bright outburst (in comparison to those flares fainter than 0.12 
crab) is about 0.65 crab estimated on August 31, 2009 (MJD 
55074), the later afterwards the higher the lower limit will be. 
The "predictions" for the 2004 and 2007 outburst (outburst 7 
and 8 in Figure [TJ based on the respective empirical relations 
derived with data before 2004 and 2007 outburst are shown as 
well. Data are from HEXTE (filled circles) and BAT (squares). 
The solid lines are the predicted values of the peak fluxes. The 
dotted lines show the corresponding prediction bounds at a level 
of 95%. 



the hard X-ray peak during the decay ph a se of t he previous out- 
burst. However as discussed in lYu et alJ d2007l) . the sum of the 
decay time of the LH state in the previous outburst and the rise 
time of the LH state in the next outburst is normally about 100- 
150 days. Therefore in reality the m inimal T w for bright out- 
bursts, as defined by lYu et all d2007l) . would be 100-150 days. 
This corresponds to F p in the range of ~ 0.04 - 0.12 crab, which 
indicates a lower limit of F p for any outburst that should be con- 
sidered in the empirical relation. This might suggest that after 
an outburst, GX 339-4 can subsequently rise up to ~ 0.12 crab 
without returning quiescence. Because of their low luminosities, 
the faint flares correspond to only a small portion of the mass in 
the disk. This consistently explains why using the empirical re- 
lation without including the faint flares can we estimate the hard 
X-ray peak flux of a bright outburst - the indicator of the disk 
mass. 

In order to get an idea on how good the estimation or pre- 
diction is, we also "predict" the hard X-ray peak flux for 2004 
and 2007 outbursts with the data before the 2004 and the 2007 
outburst respectively, and then compared the "predictions" with 
the observations (Fig [3). We then studied the deviations of the 
predicted values from the actual observed peak fluxes during 
the 2004 outburst and the 2007 outburst. The deviations are - 
0.012 crab and -0.034 crab, or 3.8% and 6.4%, respectively. 
Considering that the 10-day time binning would bring uncertain- 
ties, these predictions are extraordinarily good. The prediction 
made for the next bright outburst should have a similar accuracy. 
The hard X-ray peak of the next outburst should fall on the pre- 
diction in Fig [3] with a lower limit around 0.65 crab, which is the 
predicted hard X-ray peak flux of an outburst if it happened at 
present (around MJD 55074). 



3. DISCUSSION 

We included recent hard X-ray monitoring observations of 
GX 339-4 with Swift/BAT in addition to CGRO/BATSE and 
RXTE/HEXTE observations. We have analyzed the X-ray ob- 
servat ions of GX 339-4 in the past 18 years following I Yu et al] 
(2007) and re-examined the empirical relation between the hard 
X-ray peak flux and the outburst waiting time during bright out- 
bursts found bv lYuet al] d2007l) . We found that the hard X-ray 
peak flux of the 2007 outburst follows the empirical relation de- 
termined with observations before 2007 very well. We checked 
the potential influence of faint flares on the empirical relation. 
The empirical relation was determined based on the observations 
of bright outbursts, not including those faint flares below about 
0. 12 crab. The actually minimal waiting time required for an out- 
burst to occur consistently explains that there exists a lower limit 
of peak flux for the outburst studied here. A refined relation be- 
tween the hard X-ray peak flux and the waiting time in the past 
18 years has been obtained. Based on this relation, we can es- 
timate the hard X-ray peak flux for the next bright outburst as 
soon as it starts. It has been 750 days since the end of the most 
recent bright outburst. Based on this, we predict that the hard 
X-ray peak flux should be no less than 0.65 crab. 

One may think that during different outbursts the properties 
of the accretion flow are different, such that the radiation effi- 
ciencies differ for different outbursts while the actual mass ac- 
creted are about the same. This is not the case. The correlation 
between the hard X-ray peak flux and the peak flux of the corre- 
sponding HS state is found to hold for individua l black hole bi- 
naries and neutron star low mass X-ray binaries dYuet all [2004: 
lYu & Dolencel 120071: lYu & YanL 120091) . Given that the neutron 
star has a hard surface, the observed X-ray flux from the neutron 
star system should in general reflect the instantaneous mass ac- 
cretion rate. Therefore outbursts with different flux amplitude in 
neutron star systems should correspond to different mass accre- 
tion rate. Because the black hole systems fall on the same corre- 
lation track as those, the mass accretion rates should be different 
when GX 339-4 reaches the hard X-ray peaks during outbursts 
of different amplitudes. 

The empirical relation, confirmed by the BAT observations 
of the 2007 outburst, provides strong evidence that there is a 
link between the mass in the accretion disk and the brightest 
LH state that GX 339-4 can reach. The mechanism behind this 
link is not clear. But if the mass in the accretion disk is di- 
rectly related to the production of the hard X-ray flux, then a 
major portion of the disk should be involved in generating the 
hard X-ray flux. Independent of such accretion geometry con- 
siderations, lYu & Yarn (|2009) have recently performed a com- 
prehensive study of spectral state transitions in bright Galactic 
X-ray binaries. The results have confirmed the correlation be- 
tween LH-to-HS transition luminosity and the peak l uminosity 
of the following soft state shown in previous studies dYu et all 
120041 120071: 1 Yu & Dolencel. 120071) . and provided strong evidence 
for that: a) non-stationary accretion plays a dominant role in gen- 
erating a bright LH state and b) the rate-of-increase of the mass 
accretion rate can be the dominant parameter determining spec- 
tral state transitions. The empirical relation between the LH-to- 
HS state transition luminosity and the peak luminosity of the 
following HS state and the empirical relation studied in this pa- 
per connect the mass in the accretion disk (the cause and initial 
condition) and the peak luminosity of the hard state (the result) 
to the rate-of-increase of the mass accretion rate. This then could 
be the indicator of the initial mass which influences the overall 
development of the hard state and the soft state and the transi- 
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tions between the two. The empirical relation allows us to es- 
timate the mass in the accretion disk before an outburst in the 
special source GX 339-4. 

The phenomenon reminds us of a storage mechanism 
that works behind. This may be relevant to the phenomenon 
seen in solar flares, known as ava lanche processes (e.g., 
iLu & Hamilton 1 199 lb IWfeatlandl F2000). of which the magnetic 
field plays the major role. 
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